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Abstract. We present results of the first survey of high-redshift (< z >~ 2.3 ) O vi absorption systems along parallel lines 
of sight toward two lensed QSOs. After a careful and well-defined search, we find ten intervening O vi systems - identified 
by the presence of the /1/11031, 1037 doublet lines. Hi, and in most cases Civ, Siiv, and Cm - and eight candidate systems 
for which we do not detect H i nor other metals. We assess the veracity of these systems by applying a classification scheme. 
Within the errors, all O vi systems appear at the same redshift and have similar line strengths in front of both QSO images, 
whereas in most cases Crv or Si iv show more differences across the lines of sight, either in radial velocity or line strength. 
We conclude that (1) the coherence length of O vi must be much larger than 1 h^g kpc, and (2) an important fraction of the 
C IV absorbers may not reside in the same volume as O vi. Given the inhomogeneous character of the data -different S/N ratios 
and degrees of blending-, we pay special attention to the observational errors and their impact on the above conclusions. Since 
Doppler parameters are consistent with photoionization, we propose a model in which C iv occurs in two different photoionized 
phases, one large, with characteristic sizes of a few hundred kpc and bearing O vi, and another one a factor of ten smaller and 
containing C iii. This model is able to explain the various transverse differences observed in column density and kinematics. We 

apply the model successfully to 2 kinds of absorbers, with low and high metallicity. In the low-metallicity regime, [C/H] 2, 

we find that [C/O] ~ -0.7 is required to explain the observations, which hints at late (z S 6) rather than early metal enrichment. 
In the high-metallicity regime, the observed dissociation between O vi and C iv gas might be produced by galactic outflows. 
Altogether, the relative abundances, inhomogeneous C iv and featureless O vi are consistent with gas that has been processed 
recently before the absorption occurred (thus close to star-forming regions). Finally, we discuss briefly three associated systems 
(Zabs ~ Zem) that also show O VI. 

Key words, cosmology: observations - intergalactic medium - galaxies: halos - quasars: absorption lines - quasars: individual: 
HEI 104-1805, RXJ091I.4+055I 



1. Introduction 

Although once supposed to contain primordial material unen- 
riched by the products of star formation, the Lya forest is now 
well known to frequently exhibit associated metal lines. The 
most commonly detected metal feature at high redshift is the 
Civ doublet (e.g. Cowie & Songaila 1998; ElUson et al. 2000; 
Schaye et al. 2003; Aracil et al. 2004) due to the combination of 
its elemental abundance, convenient rest wavelength and ion- 
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ization potential. Early models which used simple recipes for 
ionization corrections and assumed uniform metallicities de- 
termined that the carbon abundance was in the range of 10"^ 
to 10"^ of the solar value (e.g Haehnelt et al. 1996; Dave et 
al 1998). More sophisticated modeling has confirmed this re- 
sult, and shown that the metallicity is most dependent on the 
gas overdensity, but largely insensitive to redshift in the range 
2 < z < 4 (Schaye et al. 2003). Although the low metalUcities 
and relatively low A'hi column densities of the associated Lya 
may hint toward an inter-galactic origin for C iv absorbers, re- 
cent work by Adelberger et al. (2003, 2005) has demonstrated 
a strong connection between the high column density C iv sys- 
tems and Lyman break galaxies at z ~ 3. Indeed, Adelberger 
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et al. (2005) infer that star-forming galaxies in their sample ex- 
hibit C IV absorption out to 80 kpc, sufficient to account for the 
bulk of observed absorbers (but note that these observations 
do not rule out an inter-galactic origin of the weak Civ ab- 
sorbers detected with pixel-statistics techniques; Fieri, Schaye 
& Aguirre 2006). 

With large datasets of high S/N spectra, extensive catalogs 
of other metal species such as Siiv (e.g., Aguirre et al. 2004) 
and Ovi (e.g., Schaye et al. 2000; Carswell, Schaye & Kim 
2002; Simcoe, Sargent & Rauch 2002, 2004; Levshakov et al. 
2003) are now becoming available. Species which represent a 
range of nucleosynthetic origins and ionization states allow us 
to probe a range of environments and physical states. For exam- 
ple, O VI is a particularly advantageous tracer of metals in the 
low density Lya forest because it has the highest ionization 
fraction for A^hi~ 10'^ cm"^ (e.g. Schaye et al. 2000; Simcoe 
et al. 2004), and because oxygen is a cosmically abundant el- 
ement. These Ovi studies have inferred enrichment of 10"^- 
IO^^Zq in IGM structures close to the mean density of the uni- 
verse at z ~ 2. However, such widespread metal pollution is 
not necessarily indicative of in-situ Population III star enrich- 
ment. Indeed, a number of authors (e.g. Aguirre et al. 2001; 
Scannapieco et al 2002; Simcoe et al. 2004) have modeled 
the contributions from Pop III stars versus galactic winds and 
concluded that the latter is the dominant enrichment source. 
The connection between O vi and galaxies has been strength- 
ened further by correlations between color-selected galaxies 
(many with spectroscopic redshifts) and absorbers by Simcoe 
et al. (2006). Ionization modeling of these O vi absorbers, many 
of which exhibit a suite of other metal species such as Nv, 
Cm and Sim, led Simcoe et al. (2006) to conclude that the 
high metallicities, large coherence scales and small cloud sizes 
were consistent with material being expelled from galaxies, e.g. 
by supernova-driven winds. However, others have argued that 
early populations of massive stars contribute a significant frac- 
tion of the metals, up to 60% according to Qian & Wasserburg 
(2005). 

Further clues to the origin and environment of Ovi ab- 
sorbers can be achieved through ionization modeling which can 
constrain densities and heating mechanisms. So far, many such 
studies have been executed at low redshift. For example, mod- 
eling the population of z < 0.3 Ovi absorbers discovered by 
HST and FUSE have revealed predominantly collisionally ion- 
ized gas at r ~ 10^ K (e.g Tripp & Savage 2000; Tripp, Savage 
& Jenkins 2000; Danforth & Shull 2005; Stocke et al. 2006), 
although multi-phase models are often required (Danforth et al. 
2005). The majority of low redshift O vi absorbers have there- 
fore been associated with inter-galactic gas that is shock heated 
as it forms large scale structures, the so-called warm-hot in- 
tergalactic medium (WHIM). O vi has also now been detected 
in a large number of Galactic sightlines, both in the disk (e.g. 
Bowen et al. 2005) and in halo clouds (e.g. Savage et al. 2003), 
and in high-velocity clouds (Fox, Savage & Wakker 2006). The 
Galactic O vi absorbers have been associated with violent dy- 
namical and non-equilibrium processes (Indebetouw & Shull 
2004; Bowen et al. 2005). In conti-ast, high redshift Ovi ab- 
sorbers have so far almost universally been found to be photo- 
ionized (Bergeron et al. 2002; Carswell et al. 2002; Levshakov 



et al 2003, but see Reimers et al. 2001, and Simcoe, Sargent & 
Rauch 2002). 

The Adelberger et al. (2005) sizes of LBG metal enrich- 
ment are consistent with limits obtained from binary or grav- 
itationally lensed QSOs which can be used to determine the 
transverse size of absorbers. A number of studies (e.g., Smette 
et al. 1995; Petitjean et al. 1998; Lopez, Hagen & Reimers 
2000; Rauch, Sargent & Barlow 2001) have seen correlations 
on scales up to tens of kpc, but with fractional column density 
variation increasing up to 50% for 100 kpc separations (e.g. 
Fig. 1 1 in Elhson et al. 2004). 

In this paper we combine the utility of gravitational lens- 
ing with ionization modeling to further investigate the physical 
state and coherence scale of Ovi absorbers. This is the first 
time that a direct measure of O vi sizes from gravitational lens- 
ing has been achieved. We complement our coherence scale 
limits derived from lensing with models of one and two phase 
ionized media. The combination of constraints from these two 
directions allows us to derive approximate sizes (under certain 
assumptions of the multi-phase structure) and relevant abun- 
dances of C, N and O. 

The targets that we use for this work are HEl 104-1805 
{Zem = 2.32) and RXJ091 1.4-H0551 (z„„ = 2.80). These ob- 
jects have been studied extensively in the lensing literature and, 
in the case of HEl 104-1805, studies on intervening absorp- 
tion line systems also exist (of which the most pertinent to 
our study are: Smette et al. 1995; Lopez et al. 1999; Rauch, 
Sargent & Barlow 2001). For RXJ091 1.4-H0551, optical low- 
resolution spectra are reported in the discovery paper by Bade 
et al. 1997, and medium-resolution spectra more recently by 
Rauch et al. (2005); here we present the first high-resolution 
spectra. Both lens systems have transverse separations of ^ 3". 
At the redshifts of the intervening systems studied here, the 
line-of-sight (LOS) separation ranges between S - 0.2 and 
1.4 kpc, for a cosmology with H^) - 70 km s"' Mpc"', and 
(Om,Qa) = (0.3,0.7). 

The paper is organized as follows: after describing the 
observations in Sect. [21 in Sect. [3] we describe the selection 
method and the individual O vi absorption systems that were 
found, along with an assessment of sample completeness. In 
Sect. |4]we analyze observed properties of those systems across 
the lines of sight. In Sect. |5] we describe results and uncertain- 
ties from photoionization models. Finally, in Sect. |6] we com- 
pare observations and models to draw inferences on the origin 
of the high-redshift O vi gas. The conclusions are outlined in 
Sect. |7] and in the Appendix we describe the associated O vi 
systems. 

2. Observations and data reduction 

We observed both QSO systems with the VLT UV- Visual 
Echelle Spectrograph (UVES) in several runs from June 2001 
through February 2003. Spectra of each QSO system were ac- 
quired simultaneously by aligning the slit with the 2 lensed im- 
ages ("A" and "B"). This configuration is well defined for the 
double QSO HEl 104-1805, but RXJ091 1.4+0551 is a quadru- 
ple system and what we shall call its "A" spectrum is actually 
the sum of three images (Al, A2, and A3 in Burud et al. 1998) 
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Table 1. VLT UVES Observations of HEl 104-1805 and RXJ091 1.4+0551 



QSO 


UVES Setup 


Wavelength 


Exp. Time 


Observing Date 


Seeing 


FWHM^ 


S/N(At 


S/N(B)'' 






[nm] 


[s] 




[arcsec] 


km s"' 






HEl 104-1805 


Die- (390+580) 


328-445,476-684 


7 600 


JunOl 


0.7-0.9 


6.7 


30,90 


15,40 




Die- (437+860) 


376-500,667-1040 


40020 


Apr01-Jan02 


0.5-1.0 










Die- (346+564) 


305-390,460,660 


11400 


AprOl-JunOl 


0.5-0.8 








RXJ091 14+0551 


Die- (390+564) 


328-445,460-560,567-660 


43 200 


Dec02-Feb03 


0.7-1.1 


7.0 


40,60 


12,20 



" Resolution of final spectrum. 

* Median signal-to-noise ratio per pixel blue and redwards of Lya emission. 
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Fig. 1. Three-cr detection limits in terms of O vi column density 
for HEl 1 04- 1 805 and RXJ09 1 1 .4+055 1 . 

which are not resolved in normal seeing conditions. Extraction 
of each 2-dimensional spectra was carried out with our own 
pipeline, which is based on simultaneous fitting of the spatial 
profiles (see Lopez et al. 2005 for details). For both QSO sys- 
tems images A and B are separated by 3^ seeing units, so the 
the fitting procedure ran smoothly. Table [T] lists the observa- 
tional and spectrum-quality parameters. 

Before coadding orders, the spectra were rebinned to a 
common vacuum-heliocentric scale with pixel size of 0.043 A. 
The final spectral resolution in the coadded spectra was FWHM 
6.7 and 7.0 km s"', for HEl 104-1805 and RXJ091 1.4+0551, 
respectively. Finally, the coadded orders were normalized by 
dividing out the response function and QSO continuum simul- 
taneously. The continuum was estimated first in the echelle or- 
ders of the A spectra by fitting cubic splines through featureless 
spectral regions. The continuum in QSO B (the faintest QSO 
component) was a scaled version of that of QSO A. Minor 
adjustments to the B continuum had to be made to account 
for differences in the flux ratios on top of emission lines. In 
RXJ091 1.4+0551, we have fitted out the mini -broad absorp- 
tion system at z = 2.57 reported in Bade et al. (1997). 

3. Absorption systems 

3.1. Selection 

Since the main goal of this paper is to study LOS differences 
between O vi absorption systems and the spectra have disparate 
S/N ratios, the selection of the systems must be made as objec- 
tive and unbiased as possible. 



We have searched for O vi systems based on three indepen- 
dent methods, which we applied to each of the LOS separately: 

1. Identify all Lye lines with A^h i ~ 10'^'^ cm"^ (checking for 
LyP and Lyy to avoid confusion with strong metal lines 
in the forest), then look for O vi at similar redshift within 
±200 km s When an Ovi candidate was found, check 
whether doublet lines could be Ly^S interlopers. Note that 
this A^H I cutoff includes systems for which O vi normally 
has not been detected by previous surveys. One caveat of 
the method is that it will fail to find highly ionized systems 
with no H i. 

2. Identify C iv absorption systems and look for O vi at similar 
redshift. 

3. Blindly search for O vi features taking into account the dou- 
blet ratio and possible forest blending. This method is lim- 
ited to clean and high S/N parts of the forest, but it is 
advantageous in finding potential broad and shallow O vi 
absorbers. This search was performed down to z > 2.0 
(HEl 104-1805 A), z > 2.1 (HEl 104-1805 B), z > 2.2 
(RXJ09 11.4+0551 A), and z > 2.5 (RXJ091 1.4+0551 B), 
which correspond to 3cr detection hmits for O vi of a; 10'^^ 
cm"^ (see Fig.[T]i. 

We found twenty one candidate Ovi systems. These are 
listed in Table |2l We then classified these candidates (column 
"class" in the Table) according to the following criteria: 

- Associated system: redshift is within y < 3 000 km s"' of 
the systemic QSO redshift. 

- Class 1: in both spectra both O vi doublet lines are detected 
at the 3cr significance level or better, and show the proper 
doublet ratios. Other metal species are present and A^hi> 
10'^-^ cm-2. 

- Class 11: like Class I, but one of the O vi doublet lines is 
detected at low significance in either of the spectra, or other 
metal species are not detected. 

- Class 111: Like Class 11 but no metals are present and A^hi< 
10'^-^ cm-2. 

All these candidate systems are shown in Figures |2] to [TTI 
and lA.ll to lA.3l Associated systems are treated in the Appendix. 

The class 111 absorbers represent the less certain identifi- 
cations, all of which are shown in Fig. |2l Almost all of these 
were discovered by method 3, except the z - 2.418109 system 
toward RXJ09 1 1 .4+055 1 . This candidate does show H i within 
±200 km s \ but the 1037 hne is also identified with 1031 of 
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Table 2. Results of search for O vi systems toward HEl 104-1805 and RXJ091 1.4+0551 



QSO 


LOS A 


LOS B 




Class 


Za 


Method" 


Zb Method" 


HE1104-1805 


2.314201 


1,2 


2.314201 


1,2 


Associated 




2.299304 


1,2 


2.299304 


1,2 


Associated 




2.200800 


1,2 


2.200822 


1,2 


II 




2.158112 


1,2 


2.158139 


1,2 


II 




2.086682 


1,2 






II 








2.052440 


2 


II 




2.036100 


3 






III 




2.016083 


3 






in 




2.010473 


2 


2.010452 


2 


II 


RXJ0911.4+0551 


2.776399 


2 


2.776399 


2 


Associated 








2.716287 


3 


III 




2.632495 


1,2 


2.632495 


1,2 


I 




2.626749 


1 






II 




2.623886 


1 






II 








2.602501 


3 


ni 




2.577051 


3 


2.573192 


3 


ni 




2.517045 


1,2 


2.516961 


2 


I 




2.437271 


1,2 


2.437265 


2 


I 




2.431602 


3 






m 




2.418109 


1,2 


2.418124 


1 


in 




2.238486 


3 






III 



" Method 1: H I as a signpost; method 2: C iv as a signpost; method 3: blind O vi search. 



the Class I system at z = 2.437271, making the former identifi- 
cation uncertain. For the rest of Class III systems, and when the 
S/N is high enough, we note that an absorption feature at the 
position of O vi appears always in the two spectra. However, 
let us emphasize that in this sub-sample there is no clear, un- 
blended case of O vi-only system. Consequently, the rest of this 
paper shall deal only with the ten remaining Class I and II sys- 
tems. 

3.2. Sample completeness 

One check of sample completeness is to compare the number 
of selected systems with previous studies. To this end we count 
systems that are class I or II (since previous studies have con- 
sidered systems with Hi). For a total redshift path surveyed of 
Az - 0.649 (based on O vi detectability in the lower S/N spec- 
tra of LOS B, which is one of the selection criteria), we find 
7 systems with W,est ^ 25 mA. The estimated number density 
of Ovi systems in our sample is thus dN/dz ~ 10.8 + 4.1 at 
< z >= 2.3. This density is consistent with previous estimates 
by e.g., Carswell, Schaye & Kim (2002; dN/dz ~ 11.4 + 4.0) 
or Reimers et al. (2001), so our particular selection criteria do 
not seem to lead to a biased sample. 

3.3. Line parameters 

We used the FITLYMAN package in MIDAS to fit Voigt pro- 
file^l to class I and II systems. When not otherwise stated, the 

' VPFIT was tried independently on several systems, yielding ex- 
actly the same results. 



O VI, C IV and H i lines were fitted independently so that iden- 
tical velocity is not a pre-requisite (except in Sect. |5] when we 
test photoionization models). Allowing the components' veloc- 
ity to be a free parameter helps quantify the extent of these dif- 
ferences along and across the LOS. The A^hi was constrained 
only in a few systems. In the much more modest S/N spectra of 
the B QSOs we attempted a fit only to O vi and C iv lines. 

In Figures [3] to [TTI we show transitions for which we were 
able to obtain line constraints. The smoothed lines are the re- 
sults of the fits and their parameters are listed in Table [3] 
The vertical dashed lines in the Figures indicate in both pan- 
els the redshift of O vi along LOS A (except in Fig. |4]l. The 
transverse (proper) distance between LOS, S{z), was calcu- 
lated using Eq. (1) of Smette et al. (1992) and Eqs. (7) and 
(8) of Tzanavaris & Carswell (2003), and assuming QSO an- 
gular separations and lens redshifts (0;ziens) = (3'.'20; 0.73) 
and (3'.'06; 0.77), respectively for HEl 104-1805 (Lidman et al. 
2000) and RXJ091 1.4+0551 (Kneib, Cohen & Hjorth 2000). 

3.4. Description of Class I and II systems 

z = 2.010473, 5 = 1.8 kpc. Class II system in HE1 104-1805, 
see Fig.|3] This system was identified through the presence of 
C IV in both spectra, and absorption features at the position of 
both Ovi doublet lines. The detection of Ovi toward QSO A 
is based on the 1037 line, albeit at low significance level. At 
the position of O vi in B the spectrum is too noisy to attempt 
any line-profile fit. C iv can be fitted in both spectra and shows 
an order of magnitude variation in column density. H i is con- 
strained in A by Lyor and Ly/? . 
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Fig. 2. Class III O vi systems toward HEl 104-1805 and RXJ091 1.4+0551 (see Table|2]). The thick black line corresponds to the 
A spectra; the thin blue one to the B spectra. 



z = 2.052440,5 = 1.5 kpc, Class II system in HE1104-1805, 
see Fig.|4] For this Ovi system unfortunately only the 1037 
lines are available. In the B spectrum, this line matches C iv and 
a weak (Nm= lO'"*^) H i line. C iii is not covered. Surprisingly 
at such small LOS separations, the H i profile is quite different 
in the A spectrum: no H i is observed at the position of O vi and 
strong H i lying redwards has no counterpart in the B spectrum. 
There is a marginal (3cr in 1548) detection of C iv toward QSO 
A at a: 80 km s"' that matches the strong H i system. 

z = 2.086682, 5 = 1.3 kpc. Class II system in HE11 04-1 805, 
see Fig. [5] In this system the detection of O vi toward QSO A 
is based on both doublet lines albeit at low significance level. 
Civ is clearly detected but blueshifted by ~ 14 km s '. Hi is 
constrained in A by Lye and Ly^S but both lines are saturated 
so we only get a lower limit on A^hi- In from a non-detection 
of Civ, we get a 5cr limit logA^(Civ) < 12.0 cm"^. At the 
position of O vi the spectrum is too noisy to attempt any line- 
profile fit. 

z = 2.1581 12, S = 0.9 kpc. Class II system in HE11 04-1 805, 
see Fig.H] For this system the detection of O vi is based on the 
1037 Une in A and the 1031 one in B (although the 1031 line 
in the latter is only a 3cr detection). In spectrum A there is a 



damped Lya system at z = 1.66 that affects the continuum 
at the position of the 1031 line. C iv is clearly detected in both 
LOS. The redshifts are almost identical (Av ~ 2 km s"' accord- 
ing to the line centroids). H i is constrained in A and B only by 
saturated Lya , so the A^hi estimates have to be taken as lower 
limits. 



z = 2.200800, S = 0.2 kpc. Class II system in HEl 104-1805, 
see Fig. [7] This is a Lyman-limit system (LLS) with rich struc- 
ture in Cm and Civ (also in Siiv, not shown in the Fig.). We 
note in the A spectrum that C m and C iv seem to follow each 
other, indicating that they arise in the same gas volume. A simi- 
lar behavior is seen in B, although with some additional blend- 
ing from interlopers. Note the v ^ 30 km s"' offset between 
Hi in A and B. The detection of Ovi is based solely on the 
1037 line in A because of Lyes' forest blending. The system 
was included because in both LOS the redshift match with C iv 
is good and the small Doppler parameters indicate metal lines. 

z = 2.437271, S = 1.4 kpc. Class I system in 
RXJ091 1.4-1-0551, see Fig. |8] In this system the two Ovi 
doublet lines are clearly present in the A spectrum. Although 
blended with shallow Lya (which have been included in the 
line fits), the doublet ratios are consistent with a positive iden- 
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Fig. 3. Class II system in HEl 104-1805. Lefthand (righthand) panels show transitions observed along sightline A (B). The 
smoothed lines show Voigt profiles with line parameters listed in Table |3] In both panels the vertical dashed line indicates 
z = 2.010473, which is the redshift of O vi in the A spectrum. 



tification. There is a good match in velocity between Civ 
and O VI. C iii is unfortunately blended and therefore not con- 
strained. From non-detection of Nv we obtain the 3cr limit 
logA^(Nv) < 12.2 cm"^. For Hi, we have considered three 
components in the fit which yields reliable H i column densi- 
ties given that the non-saturated Lyy is available. The strong 
H I velocity component closest to O vi appears blue-shifted by 
~ 100 km s The component at v = km s on the other 
hand, is weak and shallow, and has the largest Doppler param- 
eter in the sample, b -1\ km s"^ 

In B, there is a clear C iv system, which appears blueshifted 
by 43 km s ' with respect to Civ (and Ovi) in A, although 
we cannot discard that undetected weak C iv lies exactly at the 
same redshift as in A. From the Lycc profile of B (Ly/? does 
not constrain H i and Lyy is too noisy) we observe that H i does 
not vary between A and B, except for the flux excess at v = 
km s ' in LOS B. The presence of O vi in B is not conclusive; 
nevertheless a 1 -component fit gives same column density as 
in A within the formal errors, and an almost identical Zj-value. 
Inclusion of Lya blends in the fit was regarded as too compli- 
cated, given the noise level. 



z = 2.516735 and z = 2.517045, 5 = 1.0 kpc. Class I sys- 
tem in RXJ091 1 .4-H0551 , see Fig.|9] This O vi candidate also 
shows the correct line-strength ratio. We fit two velocity com- 
ponents in A, and note that C iv has also two components. The 
H I was constrained by a 2-component fit to LyyS through Lye 
, with redshifts fixed at the O vi values. C iii is also detected, 
which makes this system an excellent case for ionization di- 
agnostics. From the non-detection of Si iv we get the 3cr limit 
log A^(Si iv) < 1 1.3 cm"^; from the non-detection of N iii /1989 
we get the 3cr limit logA^(Niii) < 12.4 cm"^; and from non- 
detection of N v we get the 3o- limit log A^(N v) < 12.0 cm"^. 



In B there is absorption at the expected position of O vi but 
a 2-component fit here proved less reliable due to the lower 
S/N. We fitted a 1 -component profile and obtained the same 
value within errors as the total Ovi in A. The same proce- 
dure was applied to C iv. Therefore, the /j-values displayed in 
Table[3]have no physical meaning if this system indeed is com- 
posed of 2 velocity components along LOS B too. The Hi was 
too difficult to constrain in B again because of the high noise 
level. 
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Fig. 4. Class II system in HEl 104-1805. Panels and symbols as for Fig. [3] but for z - 2.052440, which is the redshift of O vi in 
the B spectrum. Arrows mark a possible C iv system in the A spectrum. 



z = 2.623886 and z = 2.626749, S = 0.6 kpc. Class II sys- 
tems in RXJ091 1 .4-1-0551 , see Fig. [TO] For these two sys- 
tems, the predicted position of C iv lies within spectral gaps. 
We do not detect Si iv, and C iii, N v and S vi are blended. From 
non-detection of Si iv we get the 3cr limit log A^(Si iv) < 11 .3 
cm"^. This is not proof against the identifications (neither of 
these ions is expected in an O vi phase if the gas is photoion- 
ized), but it makes the O vi identification less reliable. We note 
that H I is not aligned with O vi in either system. Although at- 
tempting Voigt profile fits in B may seem daring, they again 
show that O vi cannot differ significantly from one LOS to the 
other in neither system. 



z = 2.631480 through z = 2.634119, S = 0.5 kpc. Class I 
system in RXJ091 1 .4-1-0551 , see F\g.\U\ This is a LLS with 
a total neutral hydrogen column density in A of log A^(H i) = 
16.56 cm"^. This value comes from the unsaturated Hi /1918 
line, for which we fit 4 velocity components. This system was 
evident already in the discovery spectrum of Bade et al. (1997) 
and lies v = 11 000 km s ' from the QSO systemic redshift, 
so we regard it as intervening. Although C iv falls in the gap, 
the Si IV is evident and rich in structure. Also, lower ionization 
species are detected (not shown in the Fig.) such as Cii, Niii, 
and Si m. 



We seem to have clearly detected S vi aligned with the 
strongest Si iv component. There have been few reports of in- 
tervening S VI in the literature. Savage et al. (2005) have stud- 
ied a system at z = 0.2 with A^hi ~ 15.0 cm"^, that shows O vi, 
Ne VIII, S VI, and a wide range of low-ionization species. They 
have shown that photoionization explains the low-ionization 
species, including S vi, but is not able to reproduce the O vi. 

We have fitted 9 velocity components in each spectrum. In 
A, we note the disparity between the line centroids of O vi and 
Si IV. However, the velocities of the Si iv and H i are more sim- 
ilar (although due to higher thermal broadening we do not re- 
solve as many lines in H i as in Si iv). 

The much lower S/N at O vi in spectrum B prevents a rea- 
sonable fit. To compare both LOS we first calculated the O vi 
column density in A using both line profile fitting and the 
Apparent Optical Depth Method (AODM). The total Ovi col- 
umn density from the fits was log A^(0 vi) - 14.76 + 0.07 cm"^, 
while from the AODM we obtained log A^(0 vi) - 14.74 + 0.01 
cm"^ (integrating in the range [-150, 160] km s"', and averag- 
ing the /11031 and A1Q37 values). The good match between 
these independent estimates indicates that the lines are not 
heavily saturated. We repeated the pixel sum within the same 
velocity interval along LOS B and obtained logA^(Ovi) - 
14.66 + 0.07 cm"^, indicating that the differences in Ovi col- 
umn density between A and B are negligible. Differences in 
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Fig. 5. Class II system in HEl 104-1805. Panels and symbols as for Fig.[3] but for z = 2.086682. 



kinematical structure between A and B are also much more ev- 
ident in Si IV than in O vi. 



4. Comparisons across the lines-of-sight 

4.1. Column densities 

Fig.[T2]shows a comparison between column densities in A and 
B according to the values listed in Table [3] and excluding the 
associated systems. Column densities are summed over veloc- 
ity components when the other LOS has a value labeled with 
footnote b in the Table. It can be seen that the O vi systems do 
not present significant LOS differences (the 2 upper limits -not 
displayed in the Fig.- are consistent with no differences). In 
contrast 6 out of the 8 C iv systems significantly deviate from 
the straight line of unity-slope (for the LLS at z = 2.63 note 
that we have used Siiv as a proxy for Civ). This is also evi- 
dent in the fractional changes (lower panel of the Fig.) where 
we see that while the data is consistent with a null change in 
O VI, 6 out of the 8 C iv systems in which we detect O vi show, 
on the contrary, significant changes of up to 90 %. We have to 
keep in mind that C iv column densities are measured redwards 
of the forest and at higher S/N than O vi column densities. We 
estimate that errors in O vi column densities in the noisier B 



spectrum prevent us to assess possible changes in O vi of ; 
% or lower 



40 



In order to assess the impact of blending and S/N in the 
O VI region of the spectrum compared with the C iv region, we 
have run Monte Carlo simulations. Specifically, we synthesize 
O VI and C iv lines with random redshifts and column densities 
(within the observed ranges) and merge them with the real data. 
Two samples of spectra are created; in the first, the column 
densities are identical in the two lines of sight, in the second 
we introduce a column density difference of 0.5 dex (fractional 
change of ~ 0.7). In the latter sample, the higher column den- 
sity may be present in either LOS A or B. For 300 realizations, 
we then attempt to recover the synthesized lines by using the 
same search and fitting techniques as for the real data (the only 
exception being that line widths are kept fixed at the sim- 
ulated values, which makes error bars on the column den- 
sities smaller than for the real Unes). For the sample with 
equal column densities in the two LOS, Fig. [13] shows the 
recovered column densities in the same format as in Fig.lTIl 
We find that the recovered A^o vi and A^c iv values in the syn- 
thesized systems have fractional changes < 0.2 for 90% and 
95% of the cases, respectively, i.e., consistent with equal col- 
umn densities in the two LOS. In the sample for which the 
column densities differ by 0.5 dex in the two LOS (Fig. [14) 
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Fig. 6. Class II system in HEl 104-1805. Panels and symbols as for Fig.[3] but for z = 2.1581 12. 



the recovered fractional change is on average 0.70 for O vi 
and 0.65 for Civ, with a rather small dispersion of ~ 0.10 
for both samples. The fact that column densities are recov- 
ered with high statistical significance for both mock sam- 
ples indicates that dilferences in S/N or effects of blending 
in the Lya forest should not bias our results toward en- 
hanced coherence in O vi compared with C iv. 

It is also interesting to note that two of the three 'Class I' 
systems, regarded as more certain, show respectively low and 
high fractional changes for O vi and C iv. In summary, we con- 
clude that we have detected an excess of transverse structure in 
the C IV profiles when compared with O vi, a fact that lead us 
to conclude that a significant fraction of C iv possibly does not 
trace the same volume as the O vi-bearing gas. 

At first glance, this result seems in contradiction with pre- 
vious studies of C iv in lensed QSOs where coherence has been 
found on scales considerably larger than probed here. For ex- 
ample, in the compilation of Ellison et al. (2004), there is only 
a difference of ~ 30% for LOS separated by a few kpc. The 
fractional column density difference rises to 50% for separa- 
tions ~ 100 kpc, whereas here we find 90% variation on a 1 kpc 
scale. This apparent discrepancy can be explained by the selec- 
tion criteria of absorbers for our study. Our fundamental pre- 
requisite is that an absorber exhibit O vi, whereas the majority 
of C IV absorbers in previous samples do not show this species. 
If we consider C iv-only (as opposed to O vn-C iv) systems in 



our 2 lensed QSOs, we see a much higher level of coherence, 
with fractional variations ~ 10% (Lopez et al., in preparation). 

A second implication concerns the coherence length of the 
absorbers, i.e., the scale at which significant changes in column 
density appear Systems with A^ovi~ 10'^'^ cm"^ at z = 2.5 
must be much larger than the LOS separation scales of ^ 
kpc, given their lack of structure across the LOS. Furthermore, 
in absorbers which exhibit both C iv and O vi, the coherence 
length of the former must be smaller than the latter Since po- 
tential O VI systems occurring in just one LOS are likely not 
missing from our sample - given the agreement of our O vi 
number statistics with previous estimates -, this last conclusion 
seems generic to all O vi+C iv absorbers at high redshift. 

It is interesting to compare our transverse LOS sizes with 
those inferred from the ionization models of Simcoe et al. 
(2006). For ~ 30 absorbers with multiple transitions, Simcoe 
et al. (2006) derived metallicities and cloud sizes. They con- 
cluded that the majority of Ovi absorbers have high metallic- 
ities (up to ~ 1/3Zq) and small sizes (< 1 kpc) which inhabit 
volumes of 100 - 200 kpc around high redshift, star-forming 
galaxies. Only around 25% of the Ovi absorbers modeled by 
Simcoe et al. (2006) had sizes larger than 1 kpc. However, 
Simcoe et al. (2006) also point out that cross-section arguments 
indicate much larger dimensions; combining this with the small 
cloud sizes yielded by their ionization models led them to sug- 
gest a filamentary or sheet-like geometry. In contrast to the 
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Fig. 7. Class II system in HEl 104-1805. Panels and symbols as for Fig.[3] but for z = 2.200800 and z = 2.201069. 



mostly sub-kpc sizes inferred by Simcoe et al. (2006), we find 
that all O VI absorbers in our lensed QSOs (with LOS separa- 
tions ~ 1 kpc) appear in both LOS and with very similar col- 
umn densities. This either indicates Ovi clouds with sizes > 
1 kpc or a high filling factor of small clouds. The latter sce- 
nario was proposed by Ellison et al. (2004) to explain the co- 
herence between Mg II systems toward the triply imaged QSO 
APM08279-H5255. Although a Mg II absorber observed in one 
LOS was frequently observed in the others, the difference in 
fractional column densities could be large, and the line profiles 
were often very different. However, this does not seem to be the 
case for the Ovi absorbers. In the cases where the S/N is high 
enough that good Voigt profiles could be obtained (e.g. Figures 
|7]and|9ll, the velocity profiles are very well-matched. We there- 
fore conclude that the O vi clouds included in this sample have 
sizes larger than ~ 1 kpc. 



4.2. Velocity 

Fig. [15] shows the velocity shear between LOS A and B as 
traced by C iv and O vi. We also show the velocity differences 
between Civ and Ovi along LOS A or LOS B as a func- 
tion of LOS separation. Hi was not included in the Fig. be- 
cause line centroids are not well constrained. The small num- 
ber of systems in B for which we were able to get a Voigt 
profile precludes a full quantitative analysis on velocity shifts. 
Nevertheless, we can point out the following facts: (1) there is 
a clear offset between O vi and C iv (or Si iv) velocity compo- 
nents along LOS A in the three systems with the highest Hi 
column densities, A^hi> lO'^ cm"^; the z = 2.086682 absorber 
and the two LLSs. These systems also show differences in C iv 
or Si IV between LOS A and B (Ava-b(C iv)), apparently in- 
dependent of LOS separation. Since LLSs are believed to be 
bound to collapsed structures, Ava-b(Civ) (though not 
necessarily Avciv-Ovi(A)^ 0) are expected on kpc scales. (2) 
The 2 systems at z = 2.052440 and z = 2.437271 (both probed 
at large separations) show the largest offsets both in C iv be- 
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Fig. 8. Class I system in RXJ09 11. 4+0551. Panels and symbols as for Fig.[3l but for z = 2.437271. 



tween A and B, and between Civ and Ovi in one of the LOS. 
Interestingly, these are the 2 systems in our sample where H i 
appears "detached" from the metals. These systems are consid- 
ered in more detail below. 

Altogether, velocity offsets reinforce the notion that the 
O VI and C iv in our sample occupy (at least partially) differ- 
ent gas volumes. 

4.3. Line widths 

Considering only the "A" spectra (with better S/N), in all 5 
intervening systems where b-valnes are relatively well con- 
strained because they have single velocity components, we find 
bovi < 12 km s"', implying temperatures T < 1.4 x 10^ K. At 
this temperature and range of A^hi, models of collisionally ion- 
ized plasma (e.g., Prochaska et al. 2004) predict much less O vi 
than what we observe here. This indicates that photoionization 
is likely to be the dominant process, as has been found for the 
majority of known O vi systems at high redshift (Bergeron et 
al. 2002; Carswell, Schaye & Kim 2002). 



An important exception may be the LLS at z = 2.633, some 
of whose velocity components have ^-values ^40-50 km s ' . 
These may be artifacts due to the more complex kinematics of 
this system (we note some degree of blending with H i interlop- 
ers), or a real hint of coUisional ionization. 

Regarding line-broadening mechanisms, previous samples 
contain a mix of systems with lines dominated by either tur- 
bulence or thermal broadening. Again, considering only the 
LOS A in our sample, out of the 5 intervening systems where 
Zj-values are well constrained, 4 show clear thermal broaden- 
ing, with b(OYi) X \b(Hi) within errors (z = 2.086682, 
z = 2.158112, z = 2.623886, and z = 2.626749), and 
one (z = 2.437271) has anomalous /j-values. The system at 
z = 2.517045 has values that suggest a mild contribution 
from turbulent-motions, but the velocity components are not 
resolved, so b is less reliable there. Therefore, none of the sys- 
tems along LOS A seems to be turbulence-dominated. 

Finally, if we compare LOS, we note that 4 out of the 5 
systems where bo vi could be measured in both spectra, present 
no transverse variations in the line widths. In contrast, the (few) 
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Fig. 9. Class I system in RXJ091 1.4+0551. Panels and symbols as for Fig.[3l but for z = 2.516735 and z = 2.517045. 



/jciv-values measured are more disparate between LOS. This is 
another hint that the C iv and O vi absorbing gas may occupy 
different locations. 

4.4. Comparison with Ovi in the Galaxy 

We can compare the results of this Section with Ovi studies 
in the Galaxy. The distribution of Ovi in the Halo is shown 
to be clumpy, with fractional changes well in excess of what 
we observe at high redshift (Savage et al. 2003); homogeneous 
O VI is only shown by very small scales (a few pc; Lehner & 
Howk, 2004). But this is not the only quantitative difference 
with our high-redshift observations. First, Savage et al. (2003) 
have indicated that C iv might be kinematically less complex 
and less disturbed than Ovi in the Halo; secondly, they find 



A^(Civ)/A^(Ovi) = 0.6+0.1 (average over - 100 LOS), in con- 
trast to all high-redshift systems studied in this paper, where we 
find A^(Civ)/A?(Ovi) ^ 0.1 (excluding the LLS). Therefore, 
the systems probed by these double LOS must be of different 
nature than the Galactic ones, a fact which should not surprise 
given the different physics involved in each of them, warm-hot 
photoionized and hot collisionally-ionized gas, respectively. 

5. Photoionization models and a case study 

For the reasons outlined in Sect. |4.3| in what follows we con- 
sider the O VI gas to be photoionized. Photoionization balance 
is governed by the ionization parameter, U, the ratio of hydro- 
gen ionizing photon density to total hydrogen density, mh- In 
principle, we could constrain U by (a) assuming an ionizing 
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Fig. 10. Class II systems in RXJ09 11. 4+0551. Panels and symbols as for Fig.|3] but for z = 2.623886 and z = 2.626749. 



background spectrum, and (b) varying the gas metallicity to 
match observed A^ovi and A^civ with model. However, this ap- 
proach still requires assumption of the [C/O] ratio and that the 
two measured species arise in the same volume. In addition, 
softer spectra are degenerate with more diffuse (larger) clouds. 
This degeneracy can be broken with input from hydrodynam- 
ical simulations which correlate «h with A^hi (with a mild de- 
pendence on temperature) at least in the low-density regime 
and when hydrostatic equilibrium applies (e.g., Schaye 2001). 
However, in general, the lack of enough ionization stages of the 
same element makes detailed photoionization modeling uncer- 
tain. 

One of the O vi systems in our sample, at z = 2.5 170, is well 
suited for testing photoionization because both C iv and C iii are 
present, so no assumptions on relative abundances of different 
elements need be made. In addition, A^hi is relatively well con- 
strained by nonsaturated Lyy and Lye . Therefore, we can use 
H I, C m and C iv as "anchors" for photoionization modeling 
and investigate whether O vi can be reproduced by the models. 

5.1. A single-phase model 

We start by asking whether photoionized gas can simultane- 
ously host C III, C IV, and O vi. To test this possibility we first 
need to re-fit the data with tied velocities. The first column in 
table |4] lists the observed column densities for a best-fit red- 



shift of z = 2.517032. Note that these differ only slightly from 
those ones in Table |3] Assuming C iii and C iv coexist in the 
same volume, we build a grid of photoionization models using 
the CLOUDY package (Ferland et al. 1993), with the geom- 
etry approximated by a parallel-slab geometry. The ionizing 
agent is described by the Haardt & Madau (1996) spectrum 
of the UV background (UVB) at z = 2.5 with normalization 
log 7912 = -21.5. The spectrum includes contributions from 
QSOs and AGNs only, whose spectra are described by power 
laws F oc V"" with spectral index a = 1.8; inclusion of star- 
forming galaxies appears to have little effect on the two car- 
bon species treated here but may have an effect on Ovi (see 
Fig. 3 in Simcoe, Sargent & Rauch 2004, and Sect. l5.4l below). 
Solar abundances are taken from Grevesse & Sauval (1998), 
except for oxygen which is taken from Allende Prieto, Lambert 
& Asplund (2001). 

Fig.[T6]shows model column densities for C iii and C iv, and 
(assuming solar [C/N,0]) also for Nv, and Ovi, as a function 
of radial size S r- S ris proportional to U through fixed 7912 and 
A^Hi- The Ciii/Civ ratio determines U - 10 " which for the 
assumed value of 7912 corresponds to S - 200 kpc ( note that 
these numbers appy to a single cloud in this system). Having 
determined U, the carbon abundance can be varied to match 
the observed C iii and C iv column densities. The best match is 
found for [C/H] - -2.20. This is somewhat over the median 
but within the spread of Lya forest metallicities around [C/H] 
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Fig. 11. Class I system in RXJ09 11. 4+0551. Panels and symbols as for Fig.|3] but for z = 2.631480, z = 2.631768, z = 2.632020, 
z = 2.632495, z = 2.632770, z = 2.633161, z = 2.633486, z = 2.633766, and z = 2.6341 19. 

= -2.8 found by Simcoe, Sargent & Rauch (2004) also from With the assumed UVB intensity, the above conditions cor- 
photoionization arguments. However, it can be seen that O vi respond to a gas density nH= 10"^ ' cm"^. The nn-A^Hirelation 
is not reproduced by the model unless oxygen is overabundant yields «h~ 10"^ cm"^ (value assumes a kinetic temperature 
with respect to carbon. T ^ 1.0 x 10^ K from the newly-fitted bcn - 12.0 km s"^). 

The CLOUDY solution would match this independent value for 
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Fig. 12. Comparison of column densities between the A and B 
lines of sight (upper panel), and fractional change between A 
and B, |A^a - A^Bl/max(A^A,A^B) for Ovi and Civ in Ovi sys- 
tems. For clarity, the two upper limits on A^(0 Vi)b have not 
been displayed. The associated systems have not been included. 



Table 4. Two-phase photoionization model for cloud at z 
2.517032 



Observed 


Model 




Low High 



[C/H] 




-2.20 


S -2.10 


logt/ 




-1.15 


-0.65 


logn [cm"'] 




-3.75 


-4.3 


Size [kpc] 




20 


S 250 


[C/O] 




01 


a -0.8 


[C/N] 






^ 0.2 


logA'(Hl) 


14.88 


14.58 


14.58 


logW(Clll) 


12.50 


12.50 


« 12.00 


logWCC IV) 


12.73 


12.43 


12.43 


log ^(0 VI) 


13.92 


12.21 


13.91 


logiV(Nv) 


< 12.00 


11.60 


<, 12.00 



" Assumed. 



Fig. 13. Same as in Fig.[T2]but for simulated data. Column den- 
sities were recovered automatically from a sample of synthetic 
O VI and C iv lines with random redshifts and column densities 
and same S/N and blending levels as in the real data (see text 
for more details). The simulated column densities are identical 
in the two lines of sight. 

5.2. A multi-phase model 

Assuming homogeneity, it is difficult to think of such large 
clouds hosting C iii-C iv gas, and yet showing transverse struc- 
ture on kpc scales, as observed in Fig. [12] An overestimate of 
the size due to an incorrect UVB normalization can be ruled 
out; uncertainties inherent to this number do exist, but these do 
not reach two orders of magnitude. A more plausible situation 
might be that C iv arises in two distinct gas phases, which in 
turn would explain the offsets observed between C iv and O vi, 
since the absorption would occur in different volumes]^ 

We investigated photoionization models composed of the 
following two media: a low-ionization phase in which most 
of the observed Cm but only a fraction of the total Civ are 
present; and a high-ionization phase in which C iv and most 
of the Ovi occurs. Using the simplest approach. Hi is dis- 
tributed so as to keep the same Civ/Hi ratio in each me- 
dia. As for the C iv (and H i) fraction, very unbalanced frac- 
tions result in conditions that either approximate to the single- 
phase model described above or produce unrealistic parame- 



a factor of 2 larger clouds, a factor of 2 more intense UVB, or 
a combination of both. 



- Note that we refer to phases of photoionized gas. As we have ar- 
gued in Sect. l4.3l we do not find evidence of coUisionally ionized gas. 
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Fig. 14. Same as in Fig.[T3]but the simulated column densities 
differ by 0.5 dex (fractional change of ~ 0.7), with the higher 
column density present in either LOS A or B. 
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and O VI along LOS A or LOS B (asterisks) as a function of 
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km s"' (from line fitting). The associated systems have not been 
included. 



(in the low-ionization phase this ratio is also consistent with the 
limit for N III -see l3.4l ). 



ters (sizes of several Mpc at [M/H] = -1 to reproduce Ovi- 
only clouds). Therefore, a "50-50" distribution of the observed 
C IV appears reasonably representative of the possible variety 
of cases. Running a new grid of CLOUDY models this time 
using the modified column densities but the same approach to 
find V and Z, we determined conditions for both phases inde- 
pendently. 

The results for the (newly-fitted) z - 2.517032 system are 
displayed in Table|4] Getting parameters for the low-ionization 
phase is straightforward because we match the C iv/C iii ratio 
with the observations without assumptions on relative abun- 
dances, as in the single-phase model. This results in a cloud of 
radial size 5^ = 20 kpc at [C/H] = -2.20. A^(0 vi) is 50 times 
lower than in the high-ionization phase (assuming solar [C/O] - 
we lack of lower ionization stages of oxygen to constrain [O/H] 
in this phase), too low to be detected in our data and even for 
overabundant oxygen it would not be detectable. To model the 
high-ionization phase a [C/O] ratio has to be assumed since 
we require that C m does not contribute significantly to the ob- 
served value. That ratio is constrained imposing that the radial 
size be smaller than the maximum size allowed by the H i line 
width, if it were purely Hubble broadening, or 5 < 250 kpc for 
this particular system. To match the data, [M/H] - -2.10 and 
[C/O] = -0.8 is required (we discuss this further in the next 
sections). In addition, our limit on N v requires [C/N] ^ 0.2 dex 



5.3. Characteristic sizes along the lines of sight 

The two-phase model explains simultaneously the observed 
column density differences in C iv, the velocity shifts between 
C IV and O vi and the homogeneity of O vi on kpc scales. In 
fact, S r ~ 20 kpc sizes, if also representative of transverse di- 
mensions, perfectly produce fractional changes in C iv of the 
order of what we observe here, as demonstrated by size esti- 
mates by statistical means from observed fractional changes in 
C IV on kpc scales (Rauch, Sargent & Barlow 2001). As for the 
high-ionization phase, O vi clouds that extend over a few hun- 
dred kpc are expected to produce featureless absorption on the 
small scales probed with the present data, and, also fitting well 
to our model, studies of wider QSO pairs have found no corre- 
lation of strong C iv systems over ~ 300 kpc scales (Becker et 
al. 2004). 

Our "50-50" model is also consistent with the good match 
between C iv and C iii line profiles in all systems where C iii is 
detected (the 2 LLSs, and z = 2.5 16 systems); those cases show 
that the fraction of C iv in the low-ionization phase cannot be 
arbitrarily low. 

According to our model, if one LOS happens to pass 
through the two phases while the other LOS does so only 
through the more extended one, then column density differ- 
ences of up to 0.3 dex should arise in C iv, C iii, and H i, while 
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Fig. 16. Single-phase model for Ovi system at z = 2.517032 
. Curves are theoretical estimates of column densities as a 
function of cloud radial dimension S r for each of the relevant 
species (labeled on top or below the curves). The model as- 
sumes a metallicity [M/H] = -2.20 and solar relative abun- 
dances. The Si IV curve lies below the y-axis range. Observed 
values and their errors all are displayed at the x-axis position 
that matches C iii and C iv. 



almost the same O vi should be measured along both LOS, a sit- 
uation that is effectively observed in our sample. Regarding H i, 
although no firm constraints could be made on column densi- 
ties, there are some indications of variations of this level in the 
systems atz = 2.158112, and atz = 2.200800. 

Finally, note that we have imposed homogeneous C iv/H i 
ratios, so our model explains C iv gradients across the LOS as 
a consequence of the smaller size of low-ionization phase only, 
not because of inhomogeneous ionization conditions (within 
this phase). On the other hand, homogeneous ionization in the 
extended O vi-phase is supported by the similar line widths ob- 
served in LOS A and B. Being thermally-dominated, these sim- 
ilar widths imply the absence of significant temperature gradi- 
ents, and thus ionization must be homogeneous (at least on the 
kpc scale and with local ionizing sources being unimportant). 

5.4. Uncertainties 

The results from previous Section suffer from some model and 
measurement uncertainties. 

The first question one may raise is whether the Ovi sys- 
tems are really photoionized. In Sect. 14. 31 we have argued that 
O VI line widths of the general population in our sample suggest 
photoionized gas. In particular for the z - 2.5 17032 system, al- 
though its O VI Doppler value is not very well constrained and 
does not rule out collisional ionization (this particular measure- 



ment is subject to systematics anyway due to the unresolved ve- 
locity components), we note that both the presence of C iii and 
the narrow C iv profile hint at photoionized gas (see Prochaska 
et al. 2004). 

A second concern is what effect errors in the continuum 
can have on the CLOUDY parameters, specially given the lim- 
ited S/N of our data in this part of the forest. To explore this 
we recomputed column densities for the z - 2.517032 sys- 
tem using flux levels that were displaced artificially by Icr. 
The new column densities differ by ~ 0.08 dex. The C iii/C iv 
ratio is therefore determined with an uncertainty of at most 
~ 0.16 dex, which translates into a ~ 0.20 dex uncertainty in 
log U. This means that a wide range of possible radial sizes, 
5^- 100-500 kpc for the high-ionization and 5,. ~ 10 - 80 
kpc for the low-ionization phase, are consistent with the data 
(of course this range assumes that Jqn is accurate to this limit 
or better). Therefore, although the derived sizes are quite sen- 
sitive to changes in column density, departures from the 50-50 
model still result in two size ranges that differ by one order of 
magnitude. 

We then consider the shape and normalization of the 
UVB. As mentioned before, inclusion of star-forming galax- 
ies (Simcoe, Sargent & Rauch 2004) does not affect the region 
encompassing the C iii and C iv ionizing thresholds, but may 
lower the UVB at the Heii edge by ~ 0.3 dex. This means 
that radial sizes (derived from the C iii/C iv ratio) remain unaf- 
fected, but our result on the C/0 ratio may be affected at this 
level. Secondly, the UVB models contain also uncertainties in 
the spectral index of the quasar spectra (assumed to be pure 
power laws). The UVB hardness affects C iv and O vi differen- 
tially (Schaye et al. 2003), thus directly affecting our results on 
the C/O ratio. For instance, a spectral index of or = 1 .5 (instead 
of the assumed a - 1.8) favors a higher prediction of O vi, and 
thus a lower C/O by ^ 0.4 dex. Finally, the UVB normalization 
has been constrained with independent methods (statistics of 
metal species and proximity effect) leading to similar results, 
and is therefore subject to less important uncertainties. 

In conclusion, besides the fact that this is just a case study 
and may not be widely applicable to absorbers in general, the 
main caveat of our model seems to be the shape of the UVB, 
affecting more the derived [C/O] than the sizes. 

6. The source of metals 

6.1. Reiative abundances 

Regardless of cloud size, important inferences on relative abun- 
dances can be made. In the single-phase model depicted in 
Fig. [T6]we see that O vi and C iv are not reproduced simultane- 
ously unless [C/O]^ -0.7, and a similar situation is found for 
the two-phase model (but let us emphasize that the predicted 
[C/O] abundance is affected by the shape of the UVB: a harder 
[flatter] spectrum, say with a - 1.5, would yield [C/O] * -0.4 
only). As for nitrogen, N v should be detected according to our 
CLOUDY models but it is not, implying [C/N]^ 0.2. 

Such abundances may conform to gas that has been en- 
riched recently before the absorption occurred - in contrast to 
the yields of an early pre-enrichment. The present abundance 
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pattern is similar to direct measurements of C/H and O/H in a 
few damped Lya systems (DLA) at similar metallicity levels 
(e.g., D'Odorico & Molaro 2004). Furthermore, low [N/O] is 
found in a fraction of low-metallicity DLA, an abundance ra- 
tio that has been interpreted as an indicator of recent chemical 
enrichment (e.g., Pettini et al. 2002; Prochaska et al. 2002). 

Our result on [C/O] is also to be confronted with measure- 
ments of carbon and oxygen in Galactic metal-poor stars. For 
example, Akerman et al. (2004) have studied the metallicity 
dependence of [C/O] and confirmed the drop of [C/O] to -0.5 
at [O/H] a! -1 ([C/O] - -0.7 if non-LTE corrections are con- 
sidered, see Fabbian et al. 2005). This is roughly the regime 
we are probing in the O vi system studied in Sect. 15.21 The ab- 
sorbing gas would be enriched by a first generation of massive 
stars, just as predicted by Galactic chemical evolution models 
after completion of the first Gyr of the Halo formation. 

Alternatively, the observed gas may be chemically young 
but have been released to the IGM long-time before the absorp- 
tion occurred (say, over 2 Gyr before, z >~ 6). Though plau- 
sible, this scenario requires another ingredient, namely that the 
gas that was further processed to reach higher C/O in that time 
has not reached yet the physical volume probed by the z = 2 
observations. 

It therefore seems that one viable way to reproduce the data 
on O VI, N v, and C iv via photoionization is with CNO abun- 
dances from late rather than early injection of metals. An ex- 
tended phase of such enriched gas might be in broad line with 
[C/Si] «i -0.8 measured in the high-z IGM and in similar ab- 
sorption systems but with a different technique (Aguirre et al. 
2004). 

6.2. Ovi outflows? 

The systems at z = 2.437271 and z - 2.052440 are perhaps 
good examples of the two-phase scenario and deserve further 
analysis. Besides their particular velocity offsets (between C iv- 
O VI and H i along LOS A, and between LOS A and B in C iv 
-see Fig.O these systems also have the highest A^(0 vi)/A^(H i) 
ratios in the sample, 1 .25 and 1 .6, respectively, so both must be 
highly ionized. Systems with these A^(0 vi)/A^(H i) ratios are 
uncommon but have been reported both at high (Reimers et al. 
2001; Bergeron & Herbert-Fort 2006) and low redshift (Tripp, 
Savage & Jenkins 2000). 

The observed velocity offsets agree well with the multi- 
phase model, with a fraction of C iv in a low-ionization phase 
that is kinematically detached from a more extended Ovi 
phase. In the z - 2.437271 system along LOS B, the fraction 
of C IV that would arise in the high-ionization phase (expected 
at the position of the dashed line) is lost in the noise. In ad- 
dition, the Hi profile in B suggests a slightly lower column 
density than in A, just as expected (Sect. I5.3l l from our two- 
phase model. A similar situation is found for the z = 2.052440 
system. 

Repeating the CLOUDY simulations for these systems is 
more uncertain because C iii is not available, and thus U can 
only be constrained if a value for [C/O] is assumed. Taking the 
z = 2.437271 system along LOS A as an example, the 50-50- 



C IV distribution yields a high-ionization phase with S ,■ - 50 
kpc for [C,0/H] = -0.75 (i.e., assuming solar [C/O]). 

This C IV gas, perturbed on kpc scales (Rauch, Sargent & 
Barlow 2001) and two orders of magnitude more metal-rich 
than the average Lya forest (Simcoe, Sargent & Rauch 2004), 
is highly suggestive of gas close to overdense regions with a 
high star-formation rate. Moreover, the transverse velocity off- 
sets and column density differences of C iv and H i also suggest 
possible gas outflows with ejection velocities of a few hundred 
km s due to violent dynamical processes driven by super- 
novae explosions. 

Highly ionized, metal-rich Ovi systems detected toward 
single QSOs have been interpreted as gas outflows before (e.g., 
Aracil et al. 2004; Simcoe et al. 2006; Fox et al. 2007). The 
tens-of-km s"' shifts we detect here in Civ might be a sig- 
nature of such possible effects across the LOS (although we 
stress that these two systems have O vi line widths consistent 
with photoionization and A^hi much lower than those probed 
by the sample of Simcoe et al. 2006). The observed profiles 
could alternatively be explained by gas associated with high- 
redshift disks or extended halos, but recall from Sect. 14.41 that 
such halos would have to be quite different in nature from the 
local one. 

Another interesting point is that only Civ shows distur- 
bance on these scales, not O vi. One possible interpretation is 
that this kind of outflows do not transport highly ionized gas. 
Note that this picture does not depend on the lower significance 
of the O VI detection in the z - 2.437271 systems along LOS B, 
since O vi aligned with C iv is discarded by the data (assuming 
same O vi/C iv ratio as in A). O vi-bearing gas outflow would 
require detection of O vi aligned with C iv in B too, which is 
definitely not the case. 

6.3. Distribution of metals 

Correlations between galaxies and metal systems at z = 2 - 3 
have recently shown that O vi is always found at distances of up 
to a few hundred kpc from an observed galaxy (Adelberger et 
al. 2005; Simcoe et al. 2006). Our constraints on characteristic 
absorber sizes and relative abundances from photoionization 
models and on sizes and kinematic properties from observed 
transverse differences all fit well into a scenario of galaxy feed- 
back. The extended O vi-phase we detect here, though quite ho- 
mogeneous on kpc scales, apparently is dominated by metals 
that have been put in place only recently by Type II SNe, in 
contrast to - or on top of - an early wide spread by Pop III 
stars. 

Other observations (e.g., Adelberger et al. 2005) have indi- 
cated that gas is able to escape galaxies and pollute the IGM 
up to large distances, and this process may result in an inhomo- 
geneous metal enrichment as suggested by simulations (e.g., 
Theuns et al. 2002, Aguirre & Schaye 2005). We may have de- 
tected a new signature of this process in the disturbed C iv that 
is associated with the Ovi systems. However, the dissociation 
between C iv and O vi observed here (both along and across the 
LOS) is somewhat uncomfortable within that simple picture 
and may also indicate that the possible outflows required to 
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transport the metals have different time scales for the different 
species (and thus ionization stages) , with O vi having been "in 
place" well before Civ. On the other hand, we do not see the 
signature of shock-heated gas; if O vi systems like the one at 
z = 2.437271 are representative of such violent processes, then 
why does the gas appear predominantly photoionized? What 
fraction of the observed O vi indeed comes from galaxy feed- 
back - in contrast to arising in pre-enriched gas - remains to 
be investigated for example with similar observations of binary 
QSOs to probe larger scales. 

7. Summary and conclusions 

We have searched for O vi systems in spatially-resolved high- 
resolution spectra of two lensed QSOs. Applying a well- 
defined searching method independently to each spectrum we 
have detected Ovi in 10 intervening Hi systems. Their red- 
shifts range between z - 2.010 and 2.633, and their Ovi col- 
umn densities between 10'^'^ and 10'^'^ cm"^. We have used 
the information provided by differences and similarities across 
the lines of sight in the line parameters of various species to 
extract physical parameters via photoionization models. Given 
the inhomogeneous character of the data -different S/N ratios 
and degrees of blending-, we have paid special attention to the 
observational errors. We have also discussed the uncertainties 
of the photoionization models. In addition to this analysis, in 
the Appendix we present 3 associated systems that also show 
absorption by O vi. Our conclusions for the intervening systems 
are as follows: 

1 . Every O vi system detected in one LOS shows O vi also in 
the adjacent LOS. Since the observed number density of 
O VI is in agreement with previous studies and since H i in- 
terlopers are believed not to show strong variations on kpc- 
scales, potential Ovi systems in only one LOS are likely 
not missing from our sample. 

2. Within uncertainties, A^o VI shows little structure across the 
LOS, which probe the gas on transverse scales of 0.2-1.4 
/i^y kpc, and the data are consistent with a null fractional 
change on these scales. Unfortunately, the low S/N in the 
Lya forest of the B sightlines does not allow us to probe 
fractional changes lower than ~ 40% and the upper er- 
ror bars are in some cases consistent with high values. On 
the other hand, the more certain 'Class 1' systems all show 
low fractional changes, thus supporting an O vi coherence 
length in excess of 1 kpc. The situation for Crv is fully 
different. Six out of the 8 systems in which we detect C iv 
show a fractional change in this species that is large (up to 
90%) and significant. 

3. Civ and Ovi systematically show velocity offsets between 
line centroids while in the case where C iii is detected, it 
seems to follow C iv. There are also significant velocity off- 
sets in C IV between LOS A and LOS B, but none of the O vi 
systems show such offsets. Doppler widths indicate pho- 
toionized gas, and these also do not seem to vary between 
LOS. 

4. The observed properties of species in the 2 LOS support 
a scenario in which an important fraction of C iv may not 



reside in the same volume as O vi. This is consistent with 
photoionization models that predict clouds with unrealis- 
tically large sizes if both ions are taken into account. A 
two-phase model in which Civ arises both in a low and 
a high-ionization regime, with characteristic sizes of a few 
tens and a few hundreds of kpc, respectively, successfully 
explains the various LOS differences. 
5. The observed O vi puts constraints also on abundances, and 
photoionization requires for these clouds CNO abundances 
consistent with gas processed in galaxies recently before 
the absorption epoch. This implies that the absorbing gas is 
likely dominated by some kind of galactic feedback and/or 
it occurs not too far away from star-forming regions, al- 
though the observed O vi is of a very different nature than 
in the Galaxy. In particular, the systems at z = 2.052 and 
z - 2A31 are consistent with the lines of sight crossing a 
metal-rich galactic outflow. However, from its kinematics 
across the LOS, O vi does not seem to be transported in the 
expelled gas. 
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Appendix A: Associated systems 

The following three systems also have clear detections of O vi 
but their redshifts are within v < 3 000 km s ' of the systemic 
redshift of the QSO, so we regard them as "associated" sys- 
tems (e.g.. Barlow & Sargent 1997). For purposes of complete- 
ness, we briefly describe them here, but do not consider them in 
our photoionization models. Figures A. 1 to A. 3 show observed 
transitions of H i, C iv, and O vi plotted in velocity scale with 
respect to arbitrary redshifts. It is important to note that the 
line-of-sight separation, S , is computed in each case assuming 
that the redshift is purely cosmological, i.e., S is an upper limit 
if the absorption is produced in outflowing material. 

z = 2.300, S = 100 pc This system Hes v = 1 800 km s^' from 
the systemic redshift of HE1104-1805. We fitted 5 velocity 
components to the O vi systems along the A line-of-sight. The 
1037 transition appears underfitted due to - we believe - shal- 
low absorption by H i interlopers. Despite this caveat, a mis- 
aligment between O vi and C iv becomes clear in the 3 blue- 
most velocity components. Observations of a similar associ- 
ated system have led Ganguly et al. (2006) to postulate different 
gas phases with distinct degrees of ionization. When compar- 
ing with QSO B, C IV shows transverse structure while the O vi 
does not, which supports the scenario of different gas regions. 

This system is also observed in Si iv. Si m, and C iii (not 
shown in the Fig.) and is thus particularly well suited for study- 
ing physical conditions of the absorbing gas. Assuming the 
QSO ionizes each of the individual clouds, the size informa- 
tion provided by the two lines-of-sight can give important con- 
straints to the conditions and location of the absorbers within 
the QSO host galaxy. 

z - 2.315, S = 25 pc This system lies only v - 450 km s"' 
from the systemic redshift of HE 1 104-1805. All transitions ap- 
pear identical between the two LOS. In particular, the Ovi is 
weak and the doublet lines have the proper strength ratio. 

z = 2.777, S = 70 pc This system lies v = 2400 km s"' from 
the systemic redshift of RXJ091 1.4H-0551 and shows transi- 
tions of O VI, C IV, N v, and C m (the latter two are not shown 
in Fig. IA.3l l. The system is composed of two "clumps" of C iv 
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Table A.l. Associated O vi systems 



Ion 


z 


logMA;!! 


logA'(B)'' 


QSO 


Hi 


2.300 


16.36(0.23:1^ 


16.53(0.23)" 


HE 1104-1805 


Civ 




14.01(0.02) 


14.32(0.02) 




Ovi 




14.69(0.07) 


14.73(0.03) 




Hi 


2.314 


15.80(0.01) 


15.33(0.04) 


HE1104-1805 


Civ 




12.77(0.01) 


12.82(0.01) 




Ovi 




13.32(0.02) 


13.38(0.04) 




Hi 


2.776 


15.85(0.05) 




RXJ091 1.4+0551 


Civ 




14.10(0.03) 


14.20(0.10) 




Ovi 




14.77(0.02) 


14.75(0.03) 





" Total column densities. 
* From Lopez et al. (1999) 



clouds around v = and v = 150 km s Only for the blue- 
most of them can a clear identification of O vi be made. This 
clump, however, shows no H i, which appears only at the red- 
shift of the redmost components. In the Fig. we plot the two 
Hi (unsaturated) transitions used to constrain A^hi in A. These 
transitions are too noisy in B to perform a fit but comparison of 
low order (saturated) Lyman lines suggests the A^hi in B should 
be similar as in A. 

The fits attempted to the metal line profiles show clearly 
that the O vi doublet lines in the blue components do not show 
the proper strength ratio, and the same is true for N v and C iv. 
This effect, seen along both LOS, is only evident in the blue 
components; the two redmost Civ lines near v - 150 km s ' 
can perfectly be fitted. 

Anomalous ratios in associated systems have been inter- 
preted as produced by absorbing clouds that do not fully 
cover the background continuum emission region of the QSO 
(Barlow & Sargent 1997). However, incomplete coverage 
would be surprising in our case given the similar column den- 
sities in O VI and C iv along LOS A and B, separated by much 
larger distances than the size of the continuum source (the same 
applies if the excess flux is due to the unresolved lensed im- 
ages). The facts that the effect is observed along both LOS and 
that the red clump does not show anomalous ratios in neither 
LOS rule out two other possible explanations namely, emission 
by the lensing galaxy or a bad continuum placement. 




-200 200 -200 200 

Relative Velocity [km/s] 



Fig. A.l. Associated system in HE1104-1805 at z = 2.300. 
Panels and symbols as for Fig. [3] but for z - 2.298472, z = 
2.298920, z = 2.299304, z = 2.300127, and z = 2.3007 11. 
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Fig.A.2. Associated system in HE1104-1805 at z = 2.314. 
Panels and symbols as for Fig. |3] but forz = 2.314201. 
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Fig.A.3. Associated system in RXJ091 1.4+0551. Panels and 
symbols as for Fig. H but for z = 2.776014, z = 2.776399, and 
z = 2.776677. 
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Table 3. Intervening O vi systems toward HEl 104-1805 A and B, and RXJ091 1.4+0551 A and B 







LOS A 






LOS B 








Ion 


z 


log A' 


b 


z 


logN 


b 


Class 


QSO 


Hi 


2.010370 


14.39(0.01) 


41.3(0.9) 








II 


HEl 104-1805 


Civ 


2.010531 


11.84(0.07) 


5.5(1.6) 


2.010452 


12.88(0.01) 


7.6(0.3) 






Ovi 


2.010473 


13.39(0.21) 


9.8(6.1) 




<13.7 








Hi 








2.052476 


13.48(0.02) 


24.5(1.6) 


II 


HEl 104- 1805 


Civ 


2.053196 


11.77(0.09) 


5.2(2.1) 


2.052499 


12.84(0.02) 


8.9(0.5) 






Ovi 


2.052481 


14.01(0.06) 


8.3(1.7) 


2.052440 


13.77(0.16) 


7.4(4.5) 






Hi 


2.086492 


>15. 11(0.01) 


29.9(0.7) 








II 


HEl 104-1805 


Civ 


2.086522 


12.19(0.04) 


9.6(1.1) 




<12.0 








Ovi 


2.086682 


13.21(0.11) 


10.5(3.5) 




<13.5 








Hi 


2.158179 


> 14.70(0.02) 


26.7(0.4) 


2.158165 


> 14.80(0.05) 


27.9(0.6) 


II 


HEl 104- 1805 


Civ 


2.158093 


12.13(0.04) 


10.2(1.4) 


2.158087 


12.52(0.03) 


9.5(1.0) 






Ovi 


2.158112 


13.22(0.11) 


4.6(2.9) 


2.158139 


13.32(0.09) 


8.6(2.0) 






Hi 




1().()1{Q.12X_ 
13.67(0.03^. 






17.00(0.08)° 




II 


HEl 104- 1805 


Civ 








13.92(0.03)* 








Ovi 


2.200800 


13.00(0.20) 


4.8(3.3) 


2.200822 


13.34(0.15) 


7.2(3.5) 






Ovi 


2.201069 


13.57(0.08) 


16.7(4.4) 


2.201088 


13.44(0.15) 


12.6(5.9) 






Hi 


2.435220 


13.93(0.04) 


37.6(2.9) 








I 


RXJ09 11.4+0551 


Hi 


2.436189 


14.86(0.02) 


36.6(1.0) 












Hi 


2.437021 


13.64(0.07) 


71.1(5.8) 












Civ 


2.437208 


12.54(0.03) 


13.9(1.2) 


2.436724 


12.86(0.04) 


6.4(1.0) 






Ovi 


2.437271 


13.55(0.04) 


10.0(1.0) 


2.437265 


13.70(0.10) 


10.6(3.7) 






Hi 


2.516735 


14.82(0.03) 


59.8(2.2) 








I 


RXJ09 11.4+0551 


Cm 


2.516871 


12.18(0.20) 


26.7(0.0) 












Civ 


2.516871 


12.67(0.14) 


26.7(5.2) 












Ovi 


2.516735 


13.62(0.10) 


12.8(2.7) 












Hi 


2.517045 


14.88(0.03) 


31.9(1.3) 












Cm 


2.516995 


12.50(0.07) 


12.7(0.0) 




< 12.53 








Civ 


2.516996 


12.73(0.12) 


12.7(1.6) 


2.516916 


12.92(0.06)* 


17.6(2.6) 






Ovi 


2.517045 


13.91(0.04) 


16.7(1.5) 


2.516961 


14.06(0.07)* 


22.8(4.2) 






Hi 


2.624316 


13.84(0.77) 


27.1(3.5) 








II 


RXI09 11.4+0551 


Hi 


2.624381 


14.08(0.44) 


23.8(2.1) 












Ovi 


2.623886 


12.85(0.06) 


7.7(1.6) 


2.623860 


13.08(0.17) 


1.9(0.9) 






Hi 


2.626526 


14.49(0.01) 


44.3(0.7) 








II 


RXI09 11.4+0551 


Ovi 


2.626749 


13.18(0.04) 


12.4(1.5) 


2.626694 


13.18(0.14) 


9.6(4.7) 






Hi 


2.631842 


15.66(0.07) 


10.0(2.3) 








I 


RXJ09 11.4+0551 


Ovi 


2.631480 


13.66(0.06) 


17.3(2.5) 












Hi 


2.632487 


15.98(0.06) 


21.6(3.4) 












Ovi 


2.631768 


13.54(0.09) 


8.9(1.2) 












Ovi 


2.632020 


13.83(0.05) 


12.9(1.1) 












Hi 


2.633119 


16.25(0.04) 


21.5(2.4) 












Ovi 


2.632495 


13.91(0.15) 


47.3(13.6) 












Ovi 


2.632770 


14.43(0.04) 


52.5(4.0) 












Hi 


2.633998 


15.66(0.08) 


16.4(3.4) 












Ovi 


2.633161 


13.26(0.07) 


6.7(1.1) 












Svi 


2.633173 


12.74(0.06) 


8.5(1.5) 












Ovi 


2.633486 


13.03(0.16) 


7.8(2.1) 












Ovi 


2.633766 


13.40(0.05) 


10.1(1.2) 












Ovi 


2.634119 


13.18(0.05) 


12.7(1.6) 




14.66(0.07:1^ 
14.18(0.40)* 








Si IV 




13.88(0.02)* 













" From Lopez et al. (1999) 

* Column density integrated over all fitted velocity components. 
AODM value. 
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